The Ising Spin Glass in dimension five : link overlaps 
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Extensive simulations are made of the link overlap in five dimensional Ising Spin Glasses (ISGs) 
through and below the ordering transition. Moments of the mean link overlap distributions (the 
kurtosis and the skewness) show clear critical maxima at the ISG ordering temperature. These 
criteria can be used as efficient tools to identify a freezing transition quite generally and in any 
dimension. In the ISG ordered phase the mean link overlap distribution develops a strong two 
peak structure, with the link overlap spectra of individual samples becoming very heterogeneous. 
There is no tendency towards a "trivial" universal single peak distribution in the range of size and 
temperature covered by the data. 
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We have studied the spin and link overlaps (defined 
below, Eqs. Q] and [2]) in some detail for cubic Ising Spin 
Glasses (ISGs) with near neighbor interactions in di- 
mension five. The apparently perverse choice of dimen- 
sion five, just below the ISG upper critical dimension 
d U cd = 6, has a number of motivations. ISGs in this 
dimension have almost never been studied through sim- 
ulations, but there are precise and reliable estimates for 
the inverse critical temperatures from High Temperature 
Series Expansion (HTSE) calculations In terms of 
numbers of spins N the sizes L = 4, 6 and 8 used here are 
equivalent to three dimensional samples of size L = 10, 20 
and 32, so from the point of view of N the present sam- 
ples are "large". However, for the same N equilibration 
to criticality and beyond is much faster than in lower 
dimensions. 

We first show that at the ordering temperature the mo- 
ments of the mean link overlap distributions show char- 
acteristic peaks, reflecting the onset of correlated spin 
clusters. This phenomenon can provide both fundamen- 
tal information on the ordering process and an efficient 
tool for the accurate evaluation of critical temperatures, 
in any dimension and for any interaction distribution. 

Then in the ordered state, rather than measuring large 
numbers of samples and averaging over the obscrvablcs 
recorded, we take an alternative approach. We study a 
limited number of samples (32 for each size) and record 
explicitly the equilibrium spin overlap P(q) and link over- 
lap Q{qt) spectra for each individual sample over a suc- 
cession of inverse temperatures. Examples of individ- 
ual P(q) distributions for ISGs have been exhibited in 
the literature, whereas for the link overlaps only mean 
[Q(qe)(L,P)] distributions 0-0] but no individual sample 
Q(qi) spectra have been published as far as we are aware. 
It turns out that individual ISG samples in fact show 
strongly idiosyncratic link overlap spectra which provide 
insight into the controversial question of the structure of 
the SG ordered phase. 

The link overlap parameter [§] in ISG numerical sim- 
ulations is the bond analogue of the intensively studied 



spin overlap. In both cases two replicas (copies) A and 
B of the same physical system are first generated and 
equilibrated; updating is then continued and the "over- 
laps" between the two replicas are recorded over long 
time intervals. The spin overlap at any instant t corre- 
sponds to the fraction q(t) of spins in A and B having 
the same orientation (both up or both down), and the 
normalized overall distribution over time is written P(q). 
The link overlap corresponds to the fraction qe(t) of links 
(or bonds or edges) between spins which are either both 
satisfied or both dissatisfied in the two replicas; the nor- 
malized overall distribution over time is written Q(qi). 
The explicit definitions are 

1 N 

and 

= ^T,St{t)Sf{t)S^{t)Sf{t) (2) 

ij 

where N is the number of spins per sample and Ni the 
number of links; spins i and j are linked, denoted by ij. 
We will indicate means taken over time for a given sam- 
ple by (• • • ) and means over sets of samples by [■••]. The 
physical distinction between the information obtained 
from P(q) and Q(qe) is frequently illustrated in terms 
of a low temperature domain picture [§]. 
The Hamiltonian is as usual 

rl = - ^2 JijStSj (3) 
ij 

with the interactions the symmetric bimodal (±J) or 
Gaussian distributions normalized to (Jfj) = 1- We will 
quote inverse temperatures (3 = 1/T. (3 C = 0.3925(35) 
for the bimodal case and J3 C = 0.420(3) for the Gaussian 
from HTSE calculations [I] . Equilibration and measure- 
ment runs were performed by standard heat bath updat- 
ing (without parallel tempering) on randomly selected 
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sites. Each sample started at infinite temperature and 
was gradually cooled until it reached its final tempera- 
ture. For temperatures near T c each sample saw at least 
10 7 sweeps before reaching its final temperature. Then 
another 10 7 sweeps were run before any measurements 
took place. Normally there were about 10 sweeps be- 
tween measurements, though less for higher temperatures 
and more for lower temperatures. At each temperature 
and sample we collected 4 • 10 5 , 6 • 10 5 and 5 • 10 6 mea- 
surements for L — 8, 6 and 4 respectively. 

For the Gaussian ISG it has been shown Q that in 
equilibrium 



[( qi ){L,P)] = l-[{\U\){L,P)]/p 



(4) 



where [(\U\)(L, (3)} is the mean energy per bond. The 
Gaussian data presented here satisfy this equilibrium 
condition over the full temperature range used; the bi- 
modal samples equilibrated faster than the Gaussian ones 
and were equilibrated for as long times so we will con- 
sider that for present purposes effective equilibration has 
been reached. However we will note elsewhere that the 
condition Eq. ((4]) is necessary but not stringent enough 
for true equilibration. 

For the symmetric bimodal ISG there are simple rules 
on (qe). If p s (P) is the probability that a bond is satisfied, 



\U{0)\ = l-2p a (fi) 



(5) 



and a strict lower limit on (qi) (uncorrelated satisfied 
bond positions) is given by 

[(qt)(L,p)] >p 2 s +(l-p s ) 2 -2p 3 (l-p s ) = [U(L 7 (3) 2 }. (6) 

In the high temperature limit |£/|(/3) — > tanh(/3) so 

[(<U)(L,P)]/[(1-\U(L,I3)\/I3)]^3. (7) 

For a pure near neighbor ferromagnet (so with transla- 
tional invariance) [(qg)(L,P)] = [U(L,(3) 2 ] at all temper- 
atures. 

In addition to the mean [(<#)(/?)], important comple- 
mentary information can be obtained from the moments 
of the Q(qe) distributions. The mean Q-kurtosis can be 
defined by 



(qe - (qe))' 



and the mean Q-skewness by 



(qe - (qe))' 



Qs = 



(qe-(qe)f) / {(qe-(qe)f) 



3/2 



(8) 



(9) 



For the bimodal ISG the distributions are Gaussian at 
high temperatures (Qk = 3,Q S = 0), with peaks in the 
mean Qk(L,f3) and Q S (L,(3) at criticality, Figure 1. The 
amplitudes of the peaks decrease with increasing L; this 
could be called an "evanescent" critical phenomenon as 
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FIG. 1: (Color online) The mean Q-kurtosis Eqn. ® for the 
sets of bimodal samples with L = 4 (blue triangles), L = 6 
(red circles) and L = 8 (black squares). 



it will disappear in the thermodynamic limit. Allowing 
for a weak finite size correction term, the positions of the 
maxima j3 max for each set of peaks tend accurately to- 
wards the HTSE /3 C with increasing L. For the Gaussian 
ISG an analogous peak behavior is observed. 

Physically, the excess Q-kurtosis ("fat tailed" distri- 
butions) and Q-skewness near /3 C in ISGs must be re- 
lated to the build up of inhomogeneous temporary cor- 
related spin clusters around criticality. The data show 
that they appear when the thermodynamic limit correla- 
tion length ratio £ (/3) /L is greater than some value; then 
the Q-kurtosis and the Q-skewness each tend to a peak 
for fixed L as £(/?) diverges, so at f3 c . This simple argu- 
ment explains why in ISGs the peaks should be situated 
exactly at f3 c in the large L limit. Indeed we will show 
elsewhere that excess Q-kurtosis and Q-skcwncss peaks 
are not unique to ISGs; analogous behavior with peaks 
tending to precisely (3 C with increasing size is observed in 
a pure Ising ferromagnet, where the data are not subject 
to sample averaging noise [loj]- Potentially Q-kurtosis 
and Q-skewness measurements could be used as a stan- 
dard procedure for estimating critical temperatures. This 
method would not have the intrinsic disadvantages of the 
usual crossing point techniques. 

Turning to the ordered phase (3 > (3 C , a longstanding 
controversy concerns the correct physical description of 



ISGs for finite dimensions llj, LL2J . Generally the argu 



mcnts are focussed on the behavior to be expected at very 
large sizes and very low temperatures. Unfortunately in 
practice the combination of these limits renders direct ob- 
servations inaccessible to numerical simulations, so the 
conclusions must rely on extrapolations. Here our aim 
is to reach a phenomenological description of an ensem- 
ble of individual samples of limited size at temperatures 
somewhat below the critical temperature, so in a region 



FIG. 2: (Color online) The mean L — 6 link over- 
lap distributions [Q(qi)(L, f3)\ for inverse temperatures 
13 = 0.40,0.42,0.44,0.46,0.48,0.50,0.52,0.54,0.56,0.58,0.60 
(from left to right; olive, purple,navy, brown, yellow, pink, 
cyan, blue, green, red, black). 

which can be studied directly by simulations. 

In Figure 2 we show mean bimodal link overlap spectra 
{Q(qe)} averaged over the 32 sample set at L = 6 and (3 up 
to 0.60. Broadly similar double peak spectra have been 
observed in other ISG measurements in dimensions 3 and 
4 but the present mean spectra resemble even more 
closely mean field regime spectra in dimension 1 with al- 
gebraically decaying interactions (Lcuzzi et al Figure 2 
|6j). The presence of secondary peaks in mean [Q(<%)] 
spectra was explained by Bokil et al Q as being due 
to the coexistence of a monodomain configuration and a 
domain wall configuration, with the secondary peak cor- 
responding to one of the two replicas residing in one of 
these configuration and the other replica in the other. 
By analogy with results they had obtained on a Migdal- 
Kadanoff model, Bokil et al predicted that the secondary 
peak would melt rapidly into the main peak with increas- 
ing L. This is not what we observe. As in the mean field 
case the secondary peak becomes more and more cleanly 
separated as the size is increased and as the temperature 
is lowered. There appears to be no tendency towards 
"trivial" single peak \Q(qi)] structure in the range of L 
and P covered by the present data. 

It is of interest to go into more detail and to examine 
the behavior of the distributions P(q) and Q(qe) for in- 
dividual samples (see for instance [13[ for P(q)). Each 
of the 32 realization sets can be taken as an unbiased 
sampling of the entire population of possible realizations 
for a given size and temperature. Only a couple of these 
spectra can be exhibited here. In Figure 3 we show a 
scatter pattern of (qi) against (q 2 ) for the 32 sample set 
of L = 6 individual bimodal samples at the lowest tem- 
perature studied j3 = 0.60 ~ 1.5/3 C . It can be seen that 
(q 2 ) and (qg) are correlated linearly as would be expected 
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FIG. 3: (Color online) A scatter plot of the link overlap (qe) 
against the square of the spin overlap (q 2 ) for the set of 32 
bimodal samples with L = 6 at /3 — 0.60. 
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FIG. 4: (Color online) The temperature dependence of the 
link overlap (qe)(/3) for two extreme L = 6 bimodal samples 
: sample #22 (red circles) and sample #02 (black squares) 
together with the mean [(qi)(/3)] for the L = 6 set (upper pink 
curve) and the square of the energy per link [t/(/3) 2 ] (lower 
blue curve). The vertical line indicates the HTSE critical 
[3c = 0.3925. 

on the grounds of overlap equivalence [l4[ • We select two 
extreme samples from this set, #22 with the highest and 
#02 with the lowest average (qe). The temperature de- 
pendence of (qe)(/3) for these samples together with the 
mean [(qe){/3)] and [U((3) 2 ] are shown in Figure |U The 
degree of localization of the satisfied links can be consid- 
ered to be represented by the difference between (qt)((3) 
and the random link position limit \U({3) 2 ]. It can be 
seen that both samples begin to localize their links some- 
what before /3 C ; when the temperature is lowered further, 
the localization in sample #22 is always about twice as 
strong as that of sample #02. When the same procedure 
is followed for L — 4 and L = 8 or for the Gaussian 
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FIG. 5: (Color online) The spin overlap distribution spectra 
P(q)(/3) for the L = 6 samples #22 (red, right hand curve)and 
#02 (black, left hand curve) at fi = 0.60. (The spectra are 
symmetric about q = 0). 
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spectrum secondary peak position in Figure 2, so in the 
sense of Bokil et al it is " all wall" . It is samples of this 
type which are contributing the most strongly to the sec- 
ondary peak in the mean [Q(qi)] spectra. The only plau- 
sible interpretation which suggests itself for these spectra 
lacking EA self-overlap peaks would appear to be in terms 
of a complex configuration space, with many coexisting 
orthogonal configurations, in which A and B replicas will 
almost never find themselves in the same configuration. 

In conclusion, the characteristic critical link overlap 
moment behavior shown here in the particular case of 
dimension 5 should be observable, mutatis mutandis, at 
criticality in any dimension and for any interaction dis- 
tribution. Link overlaps can thus provide a novel and 
powerful tool for studying the critical behavior of ISGs 
and for making precise estimates of their critical temper- 
atures. 

Within the ISG ordered phase, the rich sample to 
sample heterogeneity as seen through the individual link 
overlap spectra as well as through spin overlap distri- 
butions can furnish a critical test of theoretical models. 
Mean distributions give only an incomplete view in the 
ordered phase. The present data in dimension five, close 
to but below the ISG ucd, show for some realizations of 
the interactions spin and link overlap spectra which cor- 
respond to a simple mono-domain description. For other 
realizations no self-overlap peaks are seen in the spec- 
tra, indicating that the configuration space is much more 
complex. 



FIG. 6: (Color online) The link overlap distribution spectra 
Q(qe)(P) for the L — 6 samples #22 (red, right hand curve) 
and #02 (black, left hand curve) at f3 — 0.60. The vertical 
blue line is the strict random link limit (qt{f5)) = [U(/3) 2 ]. 



interaction samples, very similar results are obtained. 

Finally we show, Figures 5 and 6, individual P(q) and 
Q{lt) spectra for these extreme samples at @ = 0.60. 
Both on the Replica Symmetry Breaking (RSB) [llj and 
the droplet [121 ] approaches one expects for individual 
spectra principal Edwards-Anderson (EA) self-overlap 
peaks at q = ±q(EA)) and qi = qe(EA), together with 
[RSB] or without [droplet] secondary peaks. The high ex- 
treme sample #22 spectra follow the simple droplet pat- 
tern and are " trivial" , but remarkably for the low extreme 
sample #02, where the P(q) spectrum is concentrated 
around q = 0, the Q(qi) spectrum shows essentially no 
weight at an EA peak position anywhere near that of the 
other spectrum. The entire Q(qe) spectrum for this sam- 
ple (and others like it) is concentrated around the mean 
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